Charge-transfer in time-dependent density-functional theory via spin-symmetry-breaking 
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Long-range charge-transfer excitations pose a major challenge for time-dependent density func- 
tional approximations. We show that spin-symmetry-breaking offers a simple solution for molecules 
composed of open-shell fragments, yielding accurate excitations at large separations when the accep- 
tor effectively contains one active electron. Unrestricted exact-exchange and self -interaction-corrected 
fimctionals are performed on one-dimensional models and the real LiH molecule within the pseu- 
dopotential approximation to demonstrate our results. 
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Although time-dependent density-functional theory 
(TDDFT) has had resounding success in predicting accu- 
rate excitation spectra in a wide variety of systems (HQ]/ 
difficulties still plague its application to certain areas. 
The problem of charge-transfer (CT) excitations has 
drawn especially significant attention in recent years 0- 
1^, due to its relevance for biomolecules, molecular con- 
ductance, solar cell design; these are systems for which 
TDDFT would be particularly attractive due to its fa- 
vorable system-size scaling. However, it has been chal- 
lenging to find a satisfactory universal solution to the 
CT problem: ab initio approaches based on modeling 
the exact kernel appear impractical, while practical ap- 
proaches tend to involve empirical parameters. Here 
we present a new approach to calculate CT excitations 
in TDDFT for certain cases, based on spin-symmetry- 
breaking. We show that accurate excitations are ob- 
tained when the acceptor is an effectively one-electron 
system, e.g. an element in Group 1 of the periodic ta- 
ble treated in the pseudopotential approximation, and 
justify why this is so. For large separations, the leading 
order behavior is captured solely from Kohn-Sham (KS) 
orbital energy differences. Results are given for model 
systems and for the LiH molecule, and suggest a type of 
Koopmans' concept for one-electron systems. 

The usual approximations in TDDFT notoriously un- 
derestimate CT excitations between fragments at large 
separation R. To leading order in 1/_R, the exact answer 
for the lowest CT frequency is: 



(1) 



where is the ionization energy of the donor, is 
the electron affinity of the acceptor and — 1/i? is the 
first electrostatic correction between the now charged 
species. (Atomic units are used throughout). It is 
well-understood why TDDFT severely underestimates 
CT 1^-6]: In TDDFT, the first step is to compute the 
Kohn-Sham (KS) orbital energy differences between oc- 
cupied {i) and unoccupied (a) orbitals, ojs = Ca — In 



a second step, these frequencies are corrected to the true 
excitations via the Hartree-exchange-correlation kernel, 
/Hxc["-o](r,r',w), which shifts and mixes the KS excita- 
tions within a matrix formulation. The kernel is a func- 
tional of the ground-state density no(r), with matrix ele- 
ments / d^r(fir'(j)^{Y)(j)a{v)fuyLc{r, r' ,u))(j)^, {Y')(j)a' (r'). For 
CT excitations, the vanishing spatial overlap at large 
separations between occupied donor and unoccupied 
acceptor orbitals sitting on different nuclei means that 
the TDDFT predictions for CT excitations reduce to the 
KS orbital energy difference, — when using usual 
semi-local functional approximations for /hxc- With 
approximate ground-state functionals the highest occu- 
pied molecular orbital (HOMO) energy, en, underesti- 
mates the true ionization energy, while the lowest un- 
occupied molecular orbital (LUMO), eL, lacks relaxation 
contributions to the electron affinity. The last few years 
have seen many methods to correct the underestima- 
tion of CT excitations, e.g. Refs. most modify the 
ground-state functional to correct the approximate KS 
homo's underestimation of /, and mix in some degree 
of Hartree-Fock, and most, but not all |@, |^ determine 
this mixing via at least one empirical parameter. Funda- 
mentally, staying within pure DFT, both the relaxation 
contributions to A and the —1/i? tail in Eq.[l]come from 
/hxc/ which must exponentially diverge with fragment 
separation fill. Worse, in the case of open-shell frag- 
ments, not covered by most of the recent fixes, addition- 
ally the exact /xc is strongly frequency-dependent |@]. 

The major reason for the awkward kernel structure in 
the case of open-shell fragments lies in the KS ground- 
state description: the HOMO (and LUMO) are delocal- 
ized over the whole molecule, quite distinct from the 
Heitler-London-like nature of the true wavefunction. In 
either the case of the exact or semi-local approxima- 
tions, their orbital energy difference tends to zero as the 
molecule is pulled apart, and so /hxc must be responsi- 
ble for the entire CT energy 0]. It is long-recognized 
that this static correlation error llQ] is the root of the 



problem of poor ground-state energies, studied exten- 
sively in molecules like H2 ITl]] , and that a simple way 
out is to allow the system to break spin-symmetry. An 
unrestricted calculation with an approximate functional 
run on a diatomic molecule, leads, at a critical internu- 
clear separation, to the spin-polarized solution obtain- 
ing the lowest ground-state energy. Albeit having incor- 
rect spin-symmetry, accurate ground-state energies are 
achieved essentially because the KS description is ren- 
dered to have one electron on each atom. 

Although less discussed, the same physics applies for 
heteroatomic molecules composed of open-shell frag- 
ments iT2ll , and suggests that symmetry -breaking could 
be a means to obtain its CT excitations. If the exact 
functional were used, the lowest-energy state remains 
correctly spin-unpolarized at any R, but at the cost 
of stark step and peak features in the bonding region, 
and strong-frequency-dependent structure in /hxc, dif- 
ficult for approximations to capture. If instead, correct 
spin-symmetry is imposed on any existing approximate 
density-functional, the groimd-state displays unphysi- 
cal fractional charges at large R, delocalized HOMO and 
LUMO orbitals, and again poor CT energies (6l [T2ll . 

The following examples show that remarkably accu- 
rate CT excitations can indeed be obtained from TDDFT 
via spin-symmetry breaking when the acceptor contains 
effectively one active electron; for large enough separa- 
tions, these are contained in simply the bare KS exci- 
tations. We present first one-dimensional models that 
enable us to compare with highly accurate numerically- 
exact solutions (computed using a Rimge-Kutta differ- 
ential equation solver as implemented in the octopus 
code 114]), and to then analyze and understand in de- 
tail why the CT in symmetry-broken TDDFT is accurate 
via the underlying potentials. For the DFT calculations 
we study the exact-excharige (EXX) and local spin den- 
sity approximation(LSD) lll3ll with self-interaction cor- 
rection (SIC). We shall see their correct long-range be- 
havior yields the correct i?-dependence at large separa- 
tions. We stress that this long-rangedness must be com- 
bined with symmetry-breaking in order to get accurate 
CT excitations: restricted calculations using these func- 
tional fail ifTsll . The functionals are implemented in oc- 
topus within the KLI approximation to OEP IIT6II . Af- 
ter studying the models, we then turn to the real LiH 
molecule. In our first model, the nuclear potentials are 
represented by short-range wells, at separation R: 

Woxt(a;) = -U^/ cosh^(a-) - U^/ cosh2(a- - R) , (2) 

where the strengths — 5.5 au, and = 6au. We 
place two electrons interacting soft-Coulombically via 

v,,{x,x') = l/^l + ix-x'y . (3) 

into this heteroatomic molecule. The unrestricted KS 
calculation yields a spin-unpolarized solution until a 
separation of about R ~ 2.3a.u. (coinciding with an 
avoided crossing in the potential energy surfaces IIT2II ) 



where it begins to break symmetry: the eigenvalues for 
the up and down spin in the ground-state begin to sep- 
arate and slowly approach those of the isolated wells as 
the molecule is pulled further apart and each electron 
settles in a different well. In contrast to Ref. 117] , local- 
ization of orbitals is in fact achieved within KLI: for two 
electrons KLI is equivalent to full OEP, but even for sys- 
tems of more than two electrons (not shown here) KLI 
yielded spin-polarized localized orbitals. 

Figure [T] plots the lowest orbital excitation energies of 
the model Eq. ll2]|-||3]l. The KS energy differences, es- 
pecially those of EXX, capture the exact CT excitations 
throughout with remarkable accuracy! 
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FIG. 1: Charge-transfer excitation energies for the model 
Eq. |2)-j3) using unrestricted SIC and EXX. The exact con- 
verge onto the asymptotic Eq. [T] (solid blue lines) soon after 
the symmetry-breaking point. (Color online) 

Why this is so can be seen by studying the underlying 
KS potentials. Consider first the limit R — > oo, where 
symmetry-breaking has placed, say, the spLn-up(down) 
electron in the left(right) well in the ground state. The 
left well is the donor for the following discussion. Fig. |2] 
plots the KS potential for the t-spin: 

(4) 

and its components at a separation of i? = 50a. u. 

In Eq. m v^-^^'' refers to the atomic acceptor(donor) 
potential (i.e. the right(left) well), and i'HM(r) = 
J (Pr'n{Y')/\Y — r'l is the Fiartree potential generated by 
density nir). As R oo, in the vicinity of the donor 
Vc,f = and Ux,t = ^Vnin^], so Ws,t = ''-'cxt- the system 
is essentially a one-electron system in this limit and lo- 
cal excitations of the f-electron are, correctly, just excited 
states of the donor. In the vicinity of the acceptor, where 
only a ^--electron lives, v^.f ^ in the limit R — > 00, and 

w^^(r - acceptor) = w^t + WhK] + Wc,t["t' "^il ■ (5) 

The Hartree potential generated by the ^-electron and 
a small correlation contribution (the two dips in the 
red curve) result in a net upward shift of Vcxt- the un- 
occupied ^-electron states living in the right-hand-well 
are shifted up in energy compared with those of the J,- 
electron (i.e. those of Vc^t)- An approximate affinity level 
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is thereby induced in the right-hand-well. So the bare 
KS orbital energy-differences yield accurate CT excita- 
tions: first, the HOMO for the t-electron is the low- 
est orbital in the left-hand- well, for which en = —I^, 
the exact ionization potential, due to Koopmans' the- 
orem i fTsi] , since both EXX and SIC are exact for one 
electron. Second, and more significantly, the LUMO 
approximates the affinity level of the right-hand-well, 
sensing the presence of the ^-electron, i.e. 
Hartree-correlation relaxation contributions to the affin- 
ity are already incorporated at the bare KS level. We 
call this eigenstate of Eq. |5j and the corresponding state 
of Ws4, for which entirely analogous analysis holds, the 
"induced affinity" levels of the right and left atoms re- 
spectively. In the limit of infinite separation, they con- 
verge onto the lowest state for the unoccupied spin of 
the isolated one-electron atom. A key point is that Ug.t 
and Vs^i are different: in any restricted calculation where 
these are the same, imoccupied levels are excitations (of 
the same A^-electron system), not affinities. 
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FIG. 2: The KS potential t;s,-|- (blue) and its components. (Color 
online). 

How well does the induced affinity level approxi- 
mate the true affinity generally? DFT theory tells us 
that the exact affinity of an iV-electron atom is A(N) = 
E{N)-E{N^\) = -eH(iV+l) where the middle expres- 
sion is computed from total energy differences while the 
third expression is computed from the highest occupied 
KS eigenvalue of the relaxed (A^+l)-electron system (see 
eg. 112]). Let us then consider the KS potential of a two- 
electron atom. Because the density is localized, an unre- 
stricted calculation yields a spin-unpolarized result. De- 
noting the two-electron ground-state density 7i2(r), 

f s ["2] = foxt + Wh ["2/2] + v\n'2\ (6) 

where Wcxt is the nuclear potential and we have noted 
that, for exact exchange, fx[«-2] = ~i^h[?^2]/2 = 
— Uh['t^2/2]. First, neglecting correlation, Eq. |6]is very 
close to Eq.|5]z/7i2 ~ In^, i.e. if, when a second electron 
is added to a well in which there is already one elec- 
tron, there is little density relaxation. This is the case 
in the model example above, since the dominant part of 
the energy of the electrons is from the external poten- 
tial. In such a case, Eqs.|5] and |6] then imply that the in- 
duced affinity level approximates the true affinity well. 



It will be a lower bound (i.e. |eL,^| < \^n[,^ = 2)|), be- 
cause electron repulsion leads to t;H["2/2] being a little 
weaker than \n\\ . This was borne out in all the EXX 
results of different models we considered. Correlation 
tends to raise the induced affinity, sometimes bringing 
it higher than the true affinity: certainly using LSD-SIC, 
i'c,t["t = 0,n| = n] forms a deeper negative well than 
ij[n2]. Shortly we will discuss examples in which the 
density relaxation is important, so that n2 is not very 
close to 2n^, and there the affinity level is not such a 
good approximation to the true affinity; consequently 
the charge-transfer excitations are not as accurate. 

The arguments above hold only for one-electron ac- 
ceptors: if the acceptor already has an electron of the 
transferring spin in it, then excitations of that spin are 
the usual constant-number excitations of TDDFT, not 
approximate affinity levels. The donor however may 
contain any number of electrons: similar models that 
have, for example, three electrons in one well and one 
in another again showed excellent CT excitations from 
the former to the latter, under spin-symmetry-breaking. 

We now extend the discussion following Eq. |4] to the 
case of an iV-electron donor and 1-electron acceptor, at 
finite but large separation. For ease of notation, assume 
again the acceptor carries a 4,-electron in the molecular 
ground state. First consider Eq. |4] near the A^-electron 
donor. For external (nuclear) potentials that decay 
Coulombically, the (TV-electron), Wg,^ + 0(l/i?^), 

where is the KS potential of the donor atom. At 
the acceptor, Us,f w;?"^ - 1/i? + 0(l/i?^) (noting that 

cancels ^'H[?^t] ~ /R, while I'x.-t- ~ — 1/i?). So, 
to leading order in R, for the t-electron, en = 
El = — 1 /R where is the induced affinity level of 
the acceptor in the infinite separation limit. Therefore, 

c^s = EL - EH = - eg - 1/^ , (7) 

as in Eq. [Ij with « — and A"^ approximated by 
the induced affinity level. (For short-ranged potentials 
as in Eq.|2l the arguments lead instead to en = + 1 
and El = f^)- Therefore, a long-ranged exchange- 
correlation as in EXX or SIC, once symmetry-broken, 
yields good CT excitations, from just its bare KS orbital 
energies, as a fimction of R, for large R. Note the impor- 
tance of correct asymptotics of the functional used for 
correct i?-dependence, as well as for accurate ionization 
potentials and induced affinities. 

For an accurate CT asymptote the density relaxation 
upon the addition of an electron must be small. But 
even when density relaxation effects are significant, 
symmetry-breaking can still be useful as we now ex- 
plain. In a practical sense, the infinite-separation limit 
itself is not so much a problem for TDDFT, because 
total ground-state energy differences computed from 
DFT 1201] can often yield reliable values for / and A 
in Eq. [T] Rather it is intermediate but large distances 
that are the challenge, where CT energies deviate from 
the asymptotic formula Eq. [1] Our symmetry-breaking 
approach can capture these deviations, going beyond 



3 



Eq. [T] The procedure is to compute the (symmetry- 
broken) KS HOMO and LUMO energy difference, but, 
when density-relaxation is large, shift it by 
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where / and A are computed from total groimd-state 
DFT energy differences. In this way, asymptotically, the 
curves approach Eq.[T]accurately, but for intermediate to 
large distances, contain correct physical deviations from 
Eq.[T]due to polarization. To illustrate this, we now con- 
sider soft-Coulomb nuclear wells: 



(9) 



and Fig. Intakes a/? = 0.7 and = 1.2. The more diffuse 
densities of these wells makes them more polarizable 
so deviations from Eq. [T] are more evident, as shown in 
Fig. 121 at intermediate separations, the (exact) CT ener- 
gies fall shy of the asymptotic Eq. [H shown as the black 
curve, due to the local polarization of the CT state to- 
wards the positive charge at the other nucleus. After 
applying the shift of Eq. |8] the unrestricted SIC results 
approach the exact results well, and capture this attrac- 
tive shift; this holds also for CT in the other direction 
(not shown). (The blue curve is the asymptote for the 
unrestricted SIC prediction, but also shifted according 
toEq.|8). 
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FIG. 3: CT from the left to the right well of model Eqs. l|9). 
Exact: 1° = 0.6206H and = 0.1199H. Unrestricted SIC: 
= eg = 0.6206H, = 0.1355H, and e^"* = 0.2183H. 
(Color online) 

So far, we have discussed CT excitations obtained 
from bare KS excitations alone and argued why these 
work so well, as demonstrated by the model examples. 
The second step of TDDFT is to apply /hxc to correct 
the KS excitations towards the exact ones; there are both 
"diagonal" terms which shift each KS excitation, as well 
as "off-diagonal" ones that mix them. For the systems 
so far discussed we expect both these effects are small, 
because (i) the diagonal term involves overlap of the oc- 
cupied and unoccupied orbitals in the excitation, which 
vanish exponentially at large distances, and (ii) there is 
little mixing with other excitations in the system. Mix- 
ing and shifting of KS excitations will be important at 
small and intermediate distances, leading to further de- 
viations from the asymptotic Eq. [TJ especially for real 
molecules, given their higher density of states. 



Turning now to a real molecule, LiH: Using a pseu- 
dopotential for the Li atom renders it an effectively one- 
electron atom, and our method captures CT in both 
directions. In Fig. |4] we plot the lowest potential en- 
ergy surfaces of the LiH molecule, computed with imre- 
stricted SIC, with the TrouUier Martins pseudopotential 
coded in octopus llT3l , compared to the highly-accurate 
configuration-interaction (CI) calculations of Ref. |2l|]. 
The induced affinity level of H is 0.0726H while that 
computed from ground-state DFT energy-differences is 
0.0264H, closer to the experimental (0.0277H), so we 
have applied the shift of Eq. [8]for R > 12. As accu- 
rate energies are unavailable for CT from H to Li, we do 
not show this curve here. The excellent agreement of the 
imrestricted approach with CI can be contrasted with re- 
stricted SIC calculations, whose collapse at smaller R is 
due to the near-degeneracy of the HOMO and LUMO 
mentioned earlier; the latter leads to convergence dif- 
ficulties for larger separations. The symmetry-broken 
SIC predicts the separation at which there is a crossing 
between the ionic curve and the Li(3s) curve very ac- 
curately, although it appears more as a direct crossing 
rather than an avoided one. 
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FIG. 4: Potential energy surfaces of the LiH molecule. Even 
the (shifted) bare KS energy-differences of unrestricted SIC are 
close to the CI ones of Ref. 121]. (Color online) 

In summary, we have shown that symmetry-breaking 
is a simple non-empirical way to obtain CT excitations 
from KS orbital energies alone, for acceptors that contain 
effectively one electron, and explained why. Strikingly 
good results were obtained for model systems as well as 
for real molecule LiH and further studies are underway. 
Applying the TDDFT kernel should improve the accu- 
racy at intermediate distances, capturing mixing of CT 
and local excitations, while these /hxc corrections will 
vanish asymptotically. 

Conceptually, the observation that for one-electron 
systems, the levels of the unoccupied spin approximate 
affinity levels can be interpreted in an extended Koop- 
mans' sense. Koopman's theorem in DFT states that 
/ = — Eh exactly, while generalized Koopman's theorem 
applies to Hartree-Fock where / « and A k. — ^l^, 

leading to the use of hybrid fimctionals for CT, men- 
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tioned earlier 10, H. Although the LUMO m exact DFT 
represents an excitation of the A^-electron system, rather 
than the {N + 1) -electron one, our results show that 
when iV = 1 in spin-DFT, the levels of the unoccu- 
pied spin can be interpreted in a generalized Koopmans' 
sense, as they approximate affinity levels. 
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